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We solve numerically a kinetic model of chiral polymerization in systems closed to matter and 
energy flow, paying special emphasis to its ability to amplify the small initial enantiomeric excesses 
due to the internal and unavoidable statistical fluctuations. The reaction steps are assumed to 
be reversible, implying a thermodynamic constraint among some of the rate constants. Absolute 
asymmetric synthesis is achieved in this scheme. The system can persist for long times in quasi- 
stationary chiral asymmetric states before racemizing. Strong inhibition leads to long-period chiral 
oscillations in the enantiomeric excesses of the longest homopolymer chains. We also calculate the 
entropy production a per unit volume and show that a increases to a peak value either before or in 
the vicinity of the chiral symmetry breaking transition. 

5-H ' PACS numbers: 05.40.Ca, 11.30.Qc, 87.15.B- 
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INTRODUCTION 



There is a growing consensus that the homochirality of biological compounds is a condition associated to life that 
should have emerged in the abiotic stages of evolution through processes of spontaneous mirror symmetry breaking 
(SMSB). This could have proceeded in a prebiotic stage, incorporating steps of increasing complexity thus leading to 
chemical systems and enantioselective chemical networks [lH3|. An important issue is therefore to identify processes 
of chirality amplification in chemical reactions. In this regard, a recent kinetic analysis of the Frank model in closed 
systems applied to the Soai reaction Q has taught us that in an actual chemical scenario, reaction networks that 
exhibit SMSB are very sensitive to chiral inductions owing to the presence of tiny initial enantiomeric excesses, 
as previously shown theoretically The stochastic scenario implies the creation of chirality from intrinsic chiral 
fluctuations and its later transmission and amplification. This can occur in far-from-equilibrium systems that undergo 
^ . dynamic phase transitions. 

The process must be coupled to others which preserve, extend, and transmit the chirality. Biological homochirality 
of living systems involves large macromolecules, therefore a central point is the relationship of the polymerization 
process with the emergence of chirality. This hypothesis has inspired recent activity devoted to modeling efforts 
aimed at understanding mirror symmetry breaking in polymerization of relevance to the origin of life. The models so 
proposed [fil-fldj are by and large, elaborate extensions and generalizations of Frank's original paradigmatic scheme 
15]. Heading this list, Sandars [6] introduced a detailed polymerization process plus the basic elements of enantiomeric 
cross inhibition as well as a chiral feedback mechanism in which only the largest polymers formed can enhance the 
production of the monomers from an achiral substrate. He treated basic numerical studies of symmetry breaking and 
bifurcation properties of this model for various values of the number of repeat units N. All the subsequent models 
cited here are variations on Sandars' original theme. Soon afterwards, Brandenburg and coworkers @ studied the 
stability and conservation properties of a modified Sandars' model and introduce a reduced N = 2 version including 
the effects of chiral bias. In [7[, they included spatial extent in this model to study the spread and propagation of chiral 
domains as well as the influence of a backround turbulent advection velocity field. The model of Wattis and Coveney 
@ differs from Sandars' in that they allow for polymers to grow to arbitrary lengths N and the chiral polymers of all 
lengths, from the dimer and upwards, act catalytically in the breakdown of the achiral source into chiral monomers. 
An analytic linear stability analysis of both the racemic and chiral solutions is carried out for the model's large N 
limit and various kinetic timescales are identified. The role of external white noise on Sandars-type polymerization 
networks including spatial extent has been explored by Gleiser and coworkers: the N = 2 truncated model introduced 
in 0] is subjected to external white noise in [Toj . chiral bias is considered in high intensity and long duration noise 
is considered in [l2| and in modified Sandars-type models with spatial extent and external noise are considered 
both for finite and infinite N, with an emphasis paid to the dynamics of chiral symmetry breaking. By contrast, Saito 
and Hyuga's [l4j model gives rise to homochiral states but differs markedly from Sandars' in that it does not invoke 
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the enantiomeric cross inhibition, allowing instead for reversibility in all the reaction steps. Their model requires 
open flow, which is the needed element of irreversibility. A different model which stands apart from the above group 
is that due to Plasson et al. |l6[ . They considered a recycled system based on reversible chemical reactions and open 
only to energy flow and without any (auto)catalytic reactions. A source of constant external energy -the element of 
irreversibility-is required to activate the monomers. This energy could be introduced into the system in a physical 
form, say, as high energy photons. A system of this kind, limited to dimerizations, was shown to have nonracemic 
stable final states for various ranges of the model parameter values and for total concentrations greater than a minimal 
value. 

The polymerization models referred to above are defined only for open flow systems which exchange matter and 
energy with the exterior. A constant source of achiral precursor is usually assumed. An unrealistic consequence is 
that homochiral chains can grow to infinite length. By contrast, most experimental procedures are carried out in 
closed and spatially bounded reaction domains and are initiated in far- from equilibrium states [H, [T^-Hil • It is thus 
crucial to have models compatible with these experimentally realistic boundary and initial conditions. The most 
immediate consequences are that polymer chains can grow to a finite maximum length and that the total system 
mass is constant. These two properties are of course intimately related. An original aspect of our work is that we 
consider the polymerization process in a closed system and include reversible reactions. This enables us to explore 
absolute asymmetric synthesis in thermodynamically closed systems (closed to matter flow) taking into account the 
backward reaction steps, which we call here "reversible reactions/reversible models" in spite of the fact that the 
values of the forward and reverse rate constants can be very different. As we are eloquently reminded by Mislow, 
mirror asymmetric states are in practice unavoidable on purely statistical grounds alone [251 ]. even in the absence 
of chiral physical forces. Absolute asymmetric synthesis is the ability of a system to amplify these statistical and 
minuscule enantiomeric biases up to observably large excesses. Thus a major goal of this work is assess the ability of 
such polymerization schemes to amplify these initial excesses, albeit if only temporarily. Asymmetric amplification is 
demonstrated to occur obeying microscopic reversibility in a reversible model closed to matter flow. This is of great 
practical interest as the racemization time scale can be substantially longer than that of the transformation of the 
initial reagents, depending on the strength of the mutual inhibition, or direct interaction between the enantiomers. 
And regarding chemical evolution, this obviates the need to invoke chiral physical fields and lends further support to 
the conviction that homochirality is a "stereochemical imperative" of molecular evolution [26| . 

The behavior of entropy in polymerization models is rarely discussed f2?| , and has not been addressed for mirror 
symmetry breaking in chiral polymerization. The entropy produced in a chemical reaction initiated out of equilibrium 
gives a measure of the dissipation during the approach to final equilibrium. In this paper, we calculate the rate of 
entropy production in chiral polymerization. Depending on the enantiomeric mutual inhibition, the entropy production 
undergoes a rapid increase to a peak value either before, or else in the vicinity of the chiral symmetry breaking 
transition, followed by an equally dramatic decrease. The system racemizes at time scales greater than that of 
polymerization, and is accompanied by a final decrease of entropy production to zero, indicating that the system has 
reached a final stable state (not to be confused with a stationary state, which must be accompanied by a nonzero 
constant entropy production). Computation of the average length homopolymer indicates that the final racemic state 
is dominated by the longest available chains. 



II. THE POLYMERIZATION MODEL 

The model we introduce and study here is modified and extended from that of Wattis and Coveney Q which 
is in turn, a generalization of Sandar's scheme @. Two salient differences that distinguish our model from these 
and other previous ones are that we (1) consider polymerization in closed systems [28[-, so that no matter 

flow is permitted with an external environment- and (2) we allow for reversible reactions in all the steps. A third 
difference is that we also include the formation (and dissociation) of the heterodimer. While heterodimer formation was 
originally contemplated in Q , it has been silently omitted from all the subsequent models 0413 that derive therefrom. 
Fragmentation in a Sandars type model has been considered previously, but was shown to yield a maximum average 
polymer length of only N = 3 repeat units (29[. 

We assume there is an achiral precursor S which can directly produce the chiral monomers L\ and R\ at a slow 
rate e as well as be consumed in processes in which homopolymers of all lengths catalyze the production of monomers. 
The specific reaction scheme we study here is composed of the following steps, where the e, (e-), k (&-), etc., denote 
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the forward (reverse) reaction rate constants and < / < 1 is the fidelity of the feedback mechanism: 



(i+/) 1(1-/) 

Li + P 

(1) 



S-Li, S + Q H_- +f) Lt + Q, S + P ^_ f) L x +P 



S f: R l7 S + P L^ +f) Ri + P, S + Q Ri + Q. 



Here 

N N 

Q = Y J nL n, P^Y. nR ^ ( 2 ) 

n— 1 n— 1 

represent a measure of the total concentrations of left-handed and right handed polymers. The amount of each chiral 
monomer produced is influenced by the total amount of chiral oligomer already present in the system. We allow for 
the monomers themselves to participate in these substrate reactions: hence N > n > 1, where N is the maximum 
chain length. The central top and bottom reactions in Eq. ([1} are cnantioselective, whereas those on the right hand 
side are non-enantioselective. The model therefore contains the features of limited enantioselectivity, first proposed 
as an alternative to the mutual inhibition of Frank [l| . 

An important observation is that differences in the Gibbs free energy AC? between initial and final states should 
be the same in all the reactions listed in Eq.([T]), which implies the thermodynamic constraint on the following forward 
and reverse reaction constants (see also [14j) 

The polymerization and chain end-termination reactions (see below) are not subject to a thermodynamic constraint. 

The monomers combine to form chirally pure polymer chains denoted by L n and R n , according to the isodesmic 
}31j stepwise reactions for 1 < n < N — 1 

kaa khb . . 

L n + L\ L n +i, R n + Ri jp R n +i, W 

and inhibition, or the chain end-termination reactions for N — 1 > n > 2: 

L n + R\ RL n , R n + L\ LR n . v 3 ) 

K ab ^ba 

These upper limits for n specified in Eqs. (|4l5p ensure that the maximum length for all oligomers produced (or 
consumed) by these reaction sets, both the homo- and heterochiral ones, is never greater than N. In the remainder of 
this paper we will consider here the natural and chiral symmetric rate assignments k aa = kbb, k a t = kb a and likewise 
for the inverse rates, k~ a — k^ b and fc~~. = k^ . There are then four independent rate constants. 
We include a separate reaction for the heterodimer formation and dissociation: 

Li+iZiJ H = L 1 R 1 , (6) 

where kh = {k a b + fc&a)/2 and k^ = (k~ b + fc ba )/2. Note that these latter two rate constants are automatically 
determined from the above choice and that L\R\ = R\L\ of course. This completes the specification of the model's 
reactions. 

The model is left-right symmetric, that is, possesses a discrete i?2 symmetry \'A{j . which is manifest in the elementary 
reaction steps, in the rate constants, and in the corresponding differential rate equations (see below). This symmetry 
can be broken spontaneously by the dynamical solutions of the differential rate equations. This model is thus apt for 
investigating spontaneous mirror symmetry breaking. Though not considered here, the effects of explicit chiral bias 
(e.g., that induced by external physical fields) can also be studied with this model by lifting the Z 2 degeneracy in 
the reaction rates, e.g., allowing for k aa ^ kbb, etc., leading to a maximum of eight independent rate constants for 
describing the reaction set in Eqs. (|4I5[) . 

Rate-equation theory as employed in chemical kinetics is used to describe the differential rate equations of the 
achiral source, the monomers, as well as the homo- and heterochiral oligomers belonging to this reaction network. 



The kinetic equations for the concentrations that follow from these elementary steps are as follows. We begin with 
the rate equations for the achiral precursor S and those corresponding to the two chiral monomers: 

^ = -2e[S\-k[S\(P + Q)+e-([L 1 ] + [R 1 ]) + ±k-[L 1 ]((l + f)Q+(l-f)P) 

+ lk-[Ri]((l + f)P + (l-f)Q), (7) 

d\T 1 k Ar_1 N ~ 1 

' '' J ^ e[S\ + ~[S]((l + f)Q + (1 - f)P) - fc 00 [£i](2[i 1 ] + i L n] + ha/kaa E ^ 



dt 



e _[i 1 ]_ii[L 1 ]((l + /)Q+(l-/)P) 



N-l N-l 



+ Ka (2[ia] + E t L «+i] + E + *fc ff > ( g ) 

= e[5] + + /)P + (1 - /)Q) - fc»[i2i](2[JJi] + E [■««] + W*» E ^ 



ra=2 n=2 

N-l N-l 

, , + + n - no^ - KJBJ^rRj + 

dt 1 1 2' 

n=2 n=l 

- e _[i2 1 ]_^l[iJ 1 ]((l + /)p + (l_/)Q) 
N-l N-l 

+ k- b (2[R 2 ] + E K+i] + W*» E [R L »]) + fc /7#- (9) 

n=2 n=2 

Whereas for iV — 1 > n > 2 we have the set of stepwise polymerization rate equations 
d[L n 



dt 
d[R n ] 



= kaa[Li]([L n -i] - [L n ]) - kablLnW + k~ b [RL n ] + fc- {[L n+l ] - [L n ]), (10) 
= fc 66 [Pi]([P„_i] - [P„]) - kbaiRn}^} + k ba [LR n ] + k b - b ([R n+1 ] - [Rn]). (11) 



dt 

Note, in accord with Eqs. (|4|5|) for the largest polymers n = N, we have instead the final pair 

d[L N " 



dt 
d[R N ] 
dt 



k aa [Li][L]\r-i] — fc~ a [Ljv], (12) 
hblRi^RN-!} - ^[Rn]. (13) 



Then the kinetic equation for the hctcrodimcr H = L\R\ (which we keep separate from the other end-chain rate 
equations): 

W. = k h [L l ][R l ]-k^[H}. (14) 

Lastly, for N — 1 > n > 2 the rate equations for the "end-spoiled" chains: 

d[LR 



dt 
d[RL ri 
dT 



= k ba [L{\[R n ] - k ba [LR n ], (15) 
= k ab [Ri][L n ] -k~ b [RL n ]. (16) 



For chemical systems closed to matter flow, the constant mass constraint that must be obeyed by the coupled system 
of differential equations at every instant is given by (the overdot denotes time derivative): 



N N-l 



[S] + 2[H] + E n([L n ] + [R n ]) + E ( n + !)([•&£»] + [LB*]) = 0. (17) 



n=1 



This relation can be verified directly using the above set of AN — 2 kinetic equations Eqs. (|7 lToT) . 
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III. ENTROPY PRODUCTION 

The entropy production rate in an irreversible process is a measure of the dissipation in that process. At equilibrium, 
the entropy production rate vanishes and is an extremum (32j . This production has been investigated recently for 
reversible versions of the Frank model [33, 34]. In those simple models, the behavior of the entropy produced near the 
chiral symmetry breaking transition as well as its subsequent temporal development, depends on whether the chemical 
system is open or closed to matter flow. We will return to this important point below. Here we consider the behavior 
of the entropy produced by polymerization reactions and monomer catalysis when the system undergoes a chiral phase 
transition as well at the later stages when the system reaches final chemical equilibrium upon racemization. 

For reactions obeying mass action kinetics, the entropy produced in any chemical reaction can be calculated straight- 
forwardly in terms of the individual elementary reaction rates [32l l35| . The rate of entropy production is the sum 
over the difference of the forward (Rj f) a nd reverse {Rjr) reaction rates multiplied by the natural logarithm of the 
ratio of the forward and reverse rates [35j : 

a(t) = R* J2(Rjf - Rjr) In (§^) > 0, (18) 

where the sum runs over each elementary reaction step j, and R* — 8.314 J mol -1 K _1 is the gas constant. Since 
our reaction scheme is set up for closed systems, equilibrium is reached after a racemization time scale t racem is 
reached, which suggests a further measure is provided by the total net entropy produced from the start of the 
reactions through chiral symmetry breaking, then on to final racemization, when the system reaches chemical and 
thermodynamic equilibrium and a(t > t racem ) = 

a T = [ a(u)du. (19) 
Jo 

This quantifies the total dissipation over the complete history of the chemical transformations under study. 

The sum in Eg. (1181) contains 2N + 4 terms. In order to determine which specific steps of the full reaction network 
provide the leading contributions to a, we group the forward and reverse reaction rates as follows: 

A. Direct monomer production 



Ri f =e[S] R lr = e4L 1 ], (20) 

R2f=e[S] R 2r = e_[R l ]. (21) 

B. Monomer catalysis 



R3f = ^(l + f)[S][Q] R 3r = ^(l + f)[L 1 ][Ql (22) 

R*t = ^ + f)[S\[P] RAr = ^(i + /)Mn (23) 

Rv = |(i - ms\[p] R 5r = ^(i - mm, (24) 

R*t = 5(1 - f)[S][Q] Rsr = y (1 - f)[R±][Q]. (25) 



For l<n<iV — lwe have 



C. Polymerization 



Rnf — k aa [Li][L n ] R% r = k- a [L n+1 ], (26) 

Rnf = hb[Rl][Rn] Rnr = Kb\Rn+l\- ( 2 7) 
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D. End-chain termination and the heterodimer 

The heterodimer rates are 

R h f = khiL^] R? = k^[H], (28) 

whereas for 2 < n < N — 1 the end-chain forward and reverse rates are given by 

R$ = k ab [L n ] [Rr] R<£ = k- b [RL n ] , (29) 

Rf f = k ba [R n ] [L x ] Rt = K a [LRn] . (30) 

IV. RESULTS 

As discussed in the Introduction, we are interested in testing the model's ability to amplify the initial small statistical 
deviations about the idealized racemic composition (25j . in systems closed to matter flow and taking microscopic 
reversibility fully into account. In order to study the sensitivity of the above reversible polymerization network to 
these minuscule initial enantiomeric excesses, a very dilute initial concentration of a scalemic (non racemic) mixture 
was employed in the calculations: the initial monomeric concentrations of [L x ]o — (1 x 10~ 6 + 1 x 10~ 15 )M and 
[i?i]o = 1 x 1CP 6 M yielding an initial chiral excess of eeo = 5 x 10~ 8 %. This is actually slightly lower [2^, [3(| 
than the excess corresponding to the initial monomer concentrations (ee = 6.1 x 10 _8 %). The initial concentration 
of the achiral substrate is [S]o — 2M, whereas those corresponding the homo- hetero-oligomers are all set to zero: 
[L n ] = [R n } = 0, for 2 < n < N, [H] = 0, and [LR n } = [RL n } = for 2 < n < N - 1. We choose e = 2x 1CT 5 , 
e_ = 10- 10 , k = 2.0, fc_ = 10- 5 , / = 0.9, k aa = k bb = 1.0, k ab - k ba = 10 3 , k~ b = fe 6 " = 1.0, fc~ = k bb = lO" 5 . For 
illustrative purposes, we consider chains that can grow to a maximum length of N = 12. This iV value is intended 
as a mean "ball-park" figure suggested by recent experiments yielding homochiral chains anywhere from N = 5 [l8T | 
up to N — 18 [24] chiral repeat units, depending on the amino acids employed and the experimental conditions. The 
differential rate equations Eqs. (|71 ll6|) were numerically integrated with the version 7 Mathematica program package 
and using a level of numerical precision, typically twenty or more significant digits, to ensure the numerical significance 
of these initial concentrations and enantiomeric excess. The results were monitored and verified to assure that total 
system mass Eq. (1171) remained constant in time. 

Results are quantified in terms of a variety of convenient chiral measures. The percent enantiomeric excess values 
of the oligomers with homochiral sequence are calculated according to (2 < n < N) 

^'Mrm* 100 ' (31) 

A global or ensemble-averaged measure of the degree of symmetry breaking is provided by the number-weighted 
enantiomeric excess rj: 

= En= 2 ([£n] - [Rn]) ^ ^ (32) 
E»= 2 ([^] + [Rn]) 

The importance of the enantiomeric excess is that it provides the order parameter for the symmetry breaking transition: 
the |ee| > is zero for chiral symmetric states and nonzero otherwise. In the latter case, the symmetry is broken. 
The average chain length of the homopolymers is given by: 



-_ En=2"([£n] + [fln]) 
El 2 ([£»] + [Rn]) ' 

and the root mean square deviation in the homochiral chain length are: 



(33) 



(n 2 ) 1/2 = ^J< (n - n) 2 > = \J < n 2 > - < n > 2 . (34) 

The angular brackets <> denote averaging with respect to the ensemble EnQ^™] + [-^n])> similar to Eq. ([33j) . It is 
important to remember that these are all time-dependent quantities. 

Temporary but rather long lived asymmetric amplification can take place, as shown in Figure[T] note the logarithmic 
time scale. The enantiomeric excess rj averaged over all chain lengths Eq. (|32[) . from the dimer on up to the maximum 
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length chain N = 12 starts off at zero value until a time on the order of t ~ 10 at which the excess increases rapidly to 
nearly 100% at SMSB. This is followed by a gradual stepwise decrease or chiral erosion characterized by the appearance 
of quasi-plateaus of approximate constancy: n falls to about 90% at t ~ 10 3 , then to about 60% at t ~ 10 6 , staying 
approximately level until the final decrease to zero occurring at a time on the order of t ~ 10 11 — 10 12 . The system 
has racemized on this time scale. No appreciable differences in 77 can be discerned when we include the monomer, that 
is, start the sums at n = 1: we still observe slow chiral erosion proceeding though quasi-steady plateaus. The rate 
of the entropy produced Eq. (|18[) exhibits an initial increase followed by a dramatic burst coinciding with the mirror 
symmetry breaking transition. This production then decreases rapidly to an exceedingly tiny but non zero value 
that remains constant during the entire period of slow chiral erosion, spanning more than ten orders of magnitude in 
time. The entropy production then goes strictly to zero when the system racemizes in complete accord with the fact 
that the system has reached chemical equilibrium. Although not displayed in Figure [TJ the substrate concentration 
falls from its initial value to zero approximately coincident with the peak structure of the entropy production, thus 
suggesting a connection between the sharp production of the latter and the change in S. The total entropy produced 
Eq. (fT9| is ut — 378.4 J mol -1 K -1 . We can identify the major contributions to the entropy production in this 
process, see Figure [2j from calculations of the partial entropies. In this way we find that the leading contribution 
comes from the monomer catalysis steps Eq.([T]), followed by the polymerization itself Eq.Q, next by the mutual 
inhibition reactions: heterodimer formation Eq.© and the "end-spoiled" cross-inhibition reactions Eq.©. The least 
important contribution comes from the direct production of monomers from the achiral substrate. The first three 
partial contributions all display a peak structure in the neighborhood of the symmetry breaking transition, with the 
corresponding peak values being displaced in time, see Figure [21 The exception to this is the entropy rate due to 
direct monomer production, which decreases monotonically. 



FIG. 1: Time evolution (logarithmic time scale) of the average enantiomeric excess n, averaged over all chains (2 < n < 12) 
Eq. (I32|l . and the associated entropy production a, Eq. (|18[) . The entropy production peaks sharply at the onset of chiral 
symmetry breaking followed by a dramatic decrease to very small values, a strictly approaches zero only at the racemization 
time scale t race m > 10 11 . Initial concentrations: [Li] = (1 x 1CT 6 + 1 x 10" 15 )M and [Ri] = 1 X 10" 6 M (ee = 5 x 10" 8 %), 
[S]o = 2M, all other homo- and heterochiral oligomer initial concentrations are zero; and for the following rates e — 2 x 10 — , 

e- = i(r 10 , k = 2.0, = i(r 5 , / = 0.9, l = k bb = 1.0, k ab = k ba = 10 3 , k~ b = k~ a = 1.0, fc- = k~ b = io~ 5 . 

A finer or more detailed measure of the degree of symmetry breaking and amplification is provided by the individual 
percent chain-length dependent enantiomeric excesses, Eq. (|31[) . A remarkable and complex dynamic behavior is 
revealed here. The time dependence of these n-dependent ee's is plotted in Figure OH note the logarithmic time scale. 
The individual ee's follow a common curve from initialization to chiral symmetry breaking, at about t ~ 10, and remain 
together at nearly 100% until about t ~ 100 at which time the common curve begins to split up into its constituents. 
Then, the percent chiral excess of each length homochiral chain behaves differently, until they again coalesce into a 
single curve upon final racemization, occurring at around t ~ 10 11 . There is a common tendency for all the ee's to 
decrease at intermediate time scales, with the largest length chains (here, N = 7, 8, 9, 10, 11, 12) passing from positive 
then to negative values of the excess. The N = 12,11,10 chains exhibit nearly —100% excess during the period 
from 10 3 to 10 4 and beyond: there has been a chiral sign reversal in the excess corresponding to the largest chains. 
This holds also for the monomer eei, plotted in the dashed curve. Except for the monomer, the individual excesses 
then all increase back to positive values at t ~ 10 6 , then from t ~ 10 7 to t ~ 10 11 , the excess increases sequentially 
as a function of the chain length n until racemization, where they all collapse to zero. The temporal behavior of 
the enantiomeric excesses of the largest chains n — 12,11,10,9,8,7 is reminiscent of strongly damped oscillations. 
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FIG. 2: Partial contributions to the total entropy production. From, (a) limited enantioselective monomer catalysis, (b) 
stepwise polymerization reactions, (c) chain-end termination and heterodimer formation reactions x 20, and (d) direct monomer 
production xlOOO. Note how the peak values of (a), (b) and (c) are displaced in time. Same initial concentrations and rate 
constants as in Fig. [1] 



This oscillatory behavior in the enantiomeric excess can be understood in terms of the evolution of the individual 
concentrations of the longest chains. To illustrate this, we focus on the time dependence of the concentrations [£12] 
and [R12], the corresponding ee experiences the largest amplitude damped oscillations, see Figure 0J For reference 
the inset diagram shows the enantiomeric excess over the entire time interval of the simulation, compare to Figure |31 
As shown in Figure HI the dominant concentration shifts back and forth between [L12] and [-R12], respectively, until 
racemization when both concentrations converge to a common value. This leads to the chiral oscillations depicted in 
the inset. The shorter chains do not experience this oscillation, as illustrated for example by the concentrations [L5] 
and [R5] plotted in Figure [SJ There, the dominant concentration is always [£5], all the way from symmetry breaking 
at t ~ 10 to racemization, at approximately t > 10 11 . The corresponding ee suffers a dip near t — 10 3 (see inset), due 
to the concentration [R5] momentarily increasing at that time, see left hand graph of Figure [5] This dip becomes more 
pronounced the longer the chain, see the sequence of curves around t ~ 10 in Figure |31 and becomes a fully-fledged 
damped chiral oscillation for the longest chains in the system. 




time(sec) 



FIG. 3: Time dependence (logarithmic scale) of the individual chain-length dependent enantiomeric excesses ee„% = + x 
100, from the start of reactions to chiral symmetry breaking, and then on to the final racemization (family of solid curves). 
Near t ~ 10 4 , the sequence of curves from top to bottom corresponds to n = 2 to n = 12, respectively. Note damped oscillatory 
behavior of the excesses corresponding to n = 7,8,9, 10, 11, 12. The dashed curve shows the chiral excess for the monomers: 
eei%. Same initial concentrations and rate constants as in Fig. [1] 



Static "snapshots" of this dynamic behavior nicely complement the evolution of the chain length dependent enan- 
tiomeric excesses. In Figure |5] we display the enantiomeric excess versus the number of chiral repeat units at selected 
time slices. In the leftmost graph, the ee's are all at 100% for all the chains. The next graph, corresponding to 
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FIG. 4: Chiral oscillations. Evolution of the individual concentrations [-L12] (the upper curve at the bifurcation) and [R12] for 
the complete time interval spanning symmetry breaking to final racemization. The inset graph gives the associated percent 
enantiomeric excess eei2 versus time. Left: expanded view of the initial stages of evolution for I <t < 10 4 . Symmetry breaking 
occurs at t ~ 10 with L12 dominating the excess. The chiral excess of these longest chains vanishes instantaneously for the first 
time at around t = 10 when the curves intersect, and then turns over such that now R12 dominates the chiral excess, see right 
hand graph (this leads to the sign flip in the excess, see inset). Right: chiral excess vanishes a second time at t ~ 10 6 when the 
two curves intersect again (compare to inset). Then from about 10 6 to 10 11 the L12 chains again dominate the chiral excess 
until racemization, when the two curves collapse to a common curve (anti-bifurcation). Same initial concentrations and rate 
constants as in Fig. [1] 
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FIG. 5: Evolution of the individual concentrations [L5] (the upper curve at the bifurcation) and [R5] for the complete time 
interval spanning symmetry breaking to final racemization. The inset graph gives the associated percent enantiomeric excess 
ees as a function of time. Left: full history of the evolution. Right: close up of the final time scales and racemization 
(anti-bifurcation) . 



t = 10 4 , shows the sign reversal tendency as a function of chain length, with the full reversal (—100%) being attained 
for the largest homochiral chains. The following graph, corresponding to t = 10 9 shows the monotone increase of 
ee with chain length. Finally, the righthand most graph shows that racemization has set in by t = 10 13 . These 
should be contrasted with Figure [31 It is interesting to point out that both the qualitative behaviors depicted in 
the second and third snapshot have been reported in two recent and independent polymerization experiments. The 
tendency of the sign reversal in ee (from positive to negative values) as a function of chain length has been observed 
in the polymerization of racemic valine (Val-NCA) and leucine (Leu-NCA)in water subject to chiral initiators (23| . 
By contrast, the monotonic increase of the percent ee with chain length has been measured in independent chiral 
amplification experiments using leucine and glycine in water [l8j starting with a 20% initial enantiomeric excess of 
the Li monomer. These static snapshots also raise the important question of when to observe the chiral amplification 
and the enantiomeric excesses. In nonlinear reaction schemes such as this one, the enantiomeric excesses one measures 
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FIG. 6: Different time slices or "snap-shots" of Fig showing dependence of the chain-length dependent enantiomeric excesses 
ee n % = x 100 for n > 2, at different selected time scales. From left to right: total chiral symmetry breaking for all 

length homopolymer chains at t = 100, next, second graph shows the sign reversal tendency for the largest chains at t = 10 4 , 
followed by the third graph, the monotone increase of chiral excess as a function if chain length at t = 10 9 , and then the fourth 
graph, the final racemization at t — 10 13 . Same initial concentrations and rate constants as in Fig. [1] 
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FIG. 7: The time evolution (logarithmic scale) of the average or mean homopolymer chain length n (upper curve) and the 
corresponding root-mean-square deviation < n 2 >5 from the mean value (lower curve). The final stable values of the mean 
and RMS values are n = 11.02 and < n 2 >5= 1.38, for t > 10 7 . This demonstrates that the final racemic mixture is dominated 
by the longer length homopolymer chains, and this is the final equilibrium configuration. Same initial concentrations and rate 
constants as in Fig. [T] 



can depend strongly on when the measurement or observation is made, that is, when one decides to terminate the 
experiment. 

Additional information regarding the homo-oligomer composition of the chemical system is provided by the average 
homochiral chain length < n >, see Eq. ([55|) . We plot this in Figure [7] along with the standard deviation about 
the mean, Eq. (|34p . The mean chain length starts off at 2, corresponding to the homodimer and then increases 
monotonically after the symmetry breaking transition, reaching a constant plateau at t = 10 7 where it remains constant 
all the way through to racemization and beyond. The final mean value < n > = 11.0. The corresponding root mean 
square indicates that the fluctuations in the mean chain length increase as the mean chain length increases but then 
drops down to a constant value (n 2 ) 1 / 2 ~ 1.4 when the average value stabilizes. This indicates that the final racemic 
composition is dominated by the longer chains: nfi na i = 11.0 ± 1.4. The racemization time scale depends on how 
"irreversible" the model is. By way of example, if we increase the rate fc_ of the reverse catalysis steps in Eq.([T]), keeping 
everything else constant, then the increased recycling of monomers back into achiral precursor S lowers this time scale 
as follows: (k-,t racem ) = (lO^ 6 , 5 x 10 12 s), (10" 5 , 5 x 10 n s), (10" 4 , 5 x 10 10 s), (10~ 3 , 1 x 10 10 s), (10~ 2 ,5 x 10 9 s). By 
the same token, if we make fc_ smaller, we can postpone racemization. 

The enantiomeric cross inhibition k a b — kt a is a determining factor in this model. By way of contrast, we consider 
a second N = 12 run with a much lower mutual inhibition than employed above, namely k a b — kb a — 20, and with the 
following inverse rates all set equal k~ a = k^ b = k~ b = k ba = 10~ 6 , but keeping the remainder of the rates as before and 
with the same initial concentrations and excess. In this situation, the symmetry breaking occurs at a later time and 
most interestingly, the entropy production now peaks well before the mirror symmetry is broken, see Figure [5] Figure 
[5] shows that the catalysis still yields the major contribution to this peak, but the second and third most important 
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FIG. 8: Time evolution (logarithmic time scale) of the average enantiomeric excess r/, averaged over all chains (2 < n < 12) 
Eq. (|32[) . and the associated entropy production a, Eq. (jTHJ) . The entropy production peaks well before (t ~ 3) the onset of 
chiral symmetry breaking (t ~ 300) followed by a dramatic decrease to very small values, a strictly approaches zero only at 
the racemization time scale t ra cem > 10 15 . Initial concentrations: [Li]o = (1 x 10" 6 + 1 x 10 _15 )M and [Ri]o = 1 x 10~ 6 M 
(eeo = 5 x 10 _8 %), [S]o = 2M, all other homo- and heterochiral oligomer initial concentrations are zero; and for the following 
rates e = 2 x 1(T 5 , e_ = 10" 10 , k = 2.0, fc_ = 10~ 5 , / = 0.9, k aa = k bb = 1.0, fe a6 = k ba = 20, Jg, = k~ a = k~ a = k M = 10" 6 




FIG. 9: Partial contributions to the total entropy production. From, (a) limited enantioselective monomer catalysis, (b) 
stepwise polymerization reactions, (c) chain-end termination and heterodimer formation reactions x 20, and (d) direct monomer 
production x 1000. Same initial concentrations and rate constants as in Fig. [8] 



contributions are now formation of end-chain spoiled oligomers followed by the polymerization, exactly opposite to 
the previous run employing the much higher mutual inhibition. The peak in a is due principally to monomer catalysis, 
and not symmetry breaking. 

The time dependence of these n-dependent ee's is plotted in Figure 1101 note the logarithmic time scale. The 
individual ee's follow a common curve from initialization to chiral symmetry breaking, at about t ~ 300, and remain 
together at nearly 100% until about t ~ 1000 at which time the common curve begins to split up into its component 
parts. Note how the shorter homopolymers tend to racemize before the longer ones, there is a sequential chiral erosion 
in the individual enantiomeric excesses that is more pronounced the shorter the chain length. This holds as well for 
the monomer, plotted in the dashed curve (contrast to the monomer behavior in Fig. [3]). Then, the percent excess of 
each length homochiral chain behaves differently, until they again coalesce into a single curve upon final racemization, 
occurring at around t ~ 10 14 . The final approach to racemization is qualitatively very similar to the case treated 
above, compare the sequence in the right hand graph of Figure [TOl to the sequence of curves in Figure [3] from roughly 
t ~ 10 7 to 10 12 . A sequence of snap-shots of the ee„'s at selected times is displayed in Figure [TT] 

Finally, we plot the average homochiral chain length < n >, see Eq. (|33|) in Figure [12] along with the standard 
deviation. The mean chain length starts off at 2, corresponding to the homodimer and then increases monotonically 
after the symmetry breaking transition, reaching a constant plateau at about t = 10 9 where it remains constant all 
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FIG. 10: Time dependence (logarithmic scale) of the individual chain-length dependent enantiomeric excesses ee n % — 
[L"] + [fl"j x fr° m the start of reactions to chiral symmetry breaking, and then on to the final racemization (family of 
solid curves). The dashed curve shows the chiral excess for the monomers. Right hand side shows a blow-up of the ee n 's for the 
time scale t = 10 9 to 10 15 , exhibiting the sequence of excesses and its final convergence to zero at racemization. Same initial 
concentrations and rate constants as in Fig. [5] 




Number of repeat units: n Number of repeat units: n Number of repeat units: n Number of repeat units: n 



FIG. 11: Different time slices or "snap-shots" of Fig[lO]showing dependence of the chain-length dependent enantiomeric excesses 
ee n % = x 100 for n > 2, at different selected time scales. From left to right: total chiral symmetry breaking for 

all length homopolymer chains at t — 10 3 , next, second graph shows the stepwise increase for the largest chains at t — 10 6 , 
followed by the third graph, the monotone increase of chiral excess as a function if chain length at t — 10 10 , and then fourth 
graph, the final racemization at approximately t = 10 14 . Same initial concentrations and rate constants as in Fig. [8] 



the way through to racemization and beyond. The final mean value < n >= 10.9. Once again, the final racemic 
composition is dominated by the longest chains. 



V. CONCLUSIONS AND DISCUSSION 



We have demonstrated that a strong chiral amplification can take place in a reversible model of chiral polymerization 
closed to matter flow and subject to constraints imposed by micro-reversibility. The inherent statistical fluctuations 
about the idealized racemic composition are modeled by an initial minuscule enantiomeric excess in a system dilute 
in the monomers. These results are important, because they suggest that spontaneous mirror symmetry breaking in 
experimental chiral polymerization can take place, and with observable and large chiral excesses without the need to 
introduce chiral initiators [24j or large initial chiral excesses (is) . Instead, the needed chiral monomers (i.e., amino 
acids) can be produced directly from an achiral precursor and amplified via catalysis. Strong mutual inhibition is 
required to amplify the initial ee to large values, very similar to what we found for the reversible Frank model in 
closed systems [4| . The chain-length dependent enantiomeric excesses depend on time in a highly nontrivial way. The 
essential rate constant is that corresponding to the enantiomeric inhibition. A most intriguing novel feature revealed 
here for appreciable enantiomeric cross inhibition is the tendency for the chain length dependent enantiomeric excesses 
to exhibit a damped oscillatory behavior before the onset of final racemization. In these conditions, the observed 
chiral excess is clearly a time dependent phenomena, though the "period" of the chiral oscillations can be quite 
long. Oscillatory dynamics in chemical reactions has been observed experimentally, and analyzed theoretically and 
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FIG. 12: The time evolution (logarithmic scale) of the average or mean homopolymer chain length n (upper curve) and the 
corresponding root-mean-square deviation < n 2 >2 from the mean value (lower curve). The final stable values of the mean 
and RMS values are n = 10.9 and < n 2 >3= 1.5, for t > 10 10 . The final racemic mixture is highly dominated by the longer 
length homopolymer chains, and this is the equilibrium configuration. Same initial concentrations and rate constants as in Fig. 

El 

numerically in simple model systems (38l . l39j ; as far as we are aware, this behavior has not been revealed previously 
as a valid dynamical solution in polymerization models. The implications for chirality transmission are far reaching: 
"memory" of the sign of the initial chiral fluctuation is washed-out by the oscillations in the enantiomeric excess, 
adding another heretofore unexpected element of randomness to the process. While the sign of the initial chiral 
fluctuation is entirely random, any subsequent chiral oscillations can further "erase" the memory of the sign of this 
initial enantiomeric excess. These oscillations cease as the system approaches its equilibrium state. Moderate values 
leads to strong temporary symmetry breaking and larger values can lead to long period damped chiral oscillations 
before final racemization takes over. 

We have also shown that the rate of entropy production per unit volume exhibits a peak value either before or 
near the vicinity of the chiral symmetry breaking transition. This increase to a peak value is mainly due to the 
catalytic production of the chiral monomers, followed next by the stepwise polymerization reactions, and then by 
the chain-end termination reactions and lastly, by the direct monomer production. The rate falls to a vanishingly 
small but constant nonzero value maintained during the intermediate time scales, then drops to zero once the system 
has racemized. Previous calculations of the entropy produced in chiral symmetry breaking transitions have been 
carried out in the Frank model. In 33], a was evaluated for a reversible open flow Frank model with a constant 
inflow of achiral substrate and a constant outflow of the mutual inhibition product. In that situation, the entropy 
production rises from small initial value and then levels off to a constant plateau after symmetry breaking. As the 
open flow keeps this system far from equilibrium, it can never racemize, and entropy is produced at a constant rate, in 
sempiternum. In [34] ]. the rate of entropy production was evaluated for a reversible open flow Frank model including 
limited enantioselectivity. Mirror symmetry is broken incompletely, and a increases sharply from the start of the 
reaction, remains at a constant level during the lifetime of the initial racemic state and then decreases to a new 
stationary value once symmetry breaking sets in. These latter authors also consider the reversible Frank model with 
constant concentration of substrate but freely varying inhibition product. The entropy production increases, a peak 
value is reached when symmetry breaking is almost complete, then decreases to very small values. On the other hand, 
for a reversible Frank model closed to matter flow, and with strong mutual inhibition, the rate of entropy production 
exhibits a sharp peak at the onset of symmetry breaking, falling to a tiny positive value until the system racemizes, 
at which time a goes to zero [37| . This latter behavior is qualitatively similar to what we find in our polymerization 
model for large inhibition. 

For sake of computational simplicity, we have considered a model wherein no generally mixed heteropolymers were 
formed, only the heterodimer LR = RL and Li? n _i and RL n ^\ for 3 < n < N. The corresponding number of 
differential equations grows linearly with the length iV as 47V — 2 . A more realistic model should of course include the 
formation of all the possible heteropolymers of a given length n, i.e., those heterochiral chains of length n containing 
r > 1 copies of L\ and s > 1 copies of Ri such that r + s = n. It is possible to build such a system, for example, 
starting with the copolymerization model of [40|, where the concentration variables are denoted as s (i)|r>i,s>o and 
c rsWlr>o,s>ii the superscript indicating the final monomer in the chain, while the double subscript r, s encodes the 



14 



number of individual L and R's making up the chain, respectively. The number of differential equations in this case 
grows quadratically with polymer length N as N(N +1), and is mathematically more involved. Detailed studies 
employing a reversible copolymerization scheme in closed systems will be presented elsewhere (4lj | . 
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